"The Magnetized Interstellar Medium" 
8-12 September 2003, Antalya, Turkey 
Eds.: B. Uyamker, W. Reich & R. Wielebinski 



Galactic Distribution and the Luminosity Function of Pulsars 

I. Yusifov and I. Kiigiik 

Department of Astronomy & Space Sciences, Faculty of Arts & Sciences, 
Erciyes University, Talas Yolu, 38039 Kayseri, Turkey 

Abstract. We have studied the radial distribution and luminosity function of normal pulsars in the Galaxy, on the 
basis of the ATNF Pulsar Catalogue where the distances are calculated according to the new electron density model 
NE200f. 

We show that the maximum of galactocentric distribution of pulsars located near 3.5 kpc and the scale-length of this 
distribution is ~ 4 kpc. The surface density of pulsars near the Galactic center (GC) is equal or slightly higher than 
that in the solar neighborhood. The local density and birth-rate of pulsars is 1.5 and 4 times higher than previous 
estimates, correspondingly before and after beaming corrections. The dependence of these results on the NE2001 model 
and recommendations for further improvement of the electron density distribution are discussed. 

We derived the luminosity function of pulsars at frequencies near 1400 MHz. In the limited band of luminosities, 
the luminosity function has a slope near —1, but a new luminosity function may be better fitted by a Log-Normal 
distribution. 

1 Introduction 

The high-frequency, sensitive Parkes Multibeam Pulsar Survey (PMPS) (Manchester et al. 2001; 
Morris et al. 2002; Kramer et al. 2003) revealed many more distant pulsars with high Dispersion 
Measures (DM). These data may allow us to investigate in more detail the statistical parameters and 
distribution of pulsars, especially in the central regions of the Galaxy, which was almost impossible in 
previous low-frequency and less-sensitive surveys. Furthermore, high-frequency searches reveal many 
pulsars at the lower as well as at the high luminosity ends of the luminosity function (LF) of pulsars. 
From the evolutionary point of view it is of great interest for the influence of such pulsars on the shape 
of pulsars' LF. 

Considering the actuality of the problem and with the appearance of publications with completely 
wrong interpretation of observational data of pulsars (Guseinov et al. 2002a, 2002b), we decided to 
estimate the density of pulsars at the GC and to take a fresh look at the radial distribution and 
luminosity function of pulsars using these new PMPS data, with the kind permission of Prof. R.N. 
Manchester, from the ATNF Pulsar Catalogue of 1412 pulsars (Manchester et al. 2002). We estimated 
the distances to the pulsars by using the NE2001 Galactic electron density model (Cordes & Lazio 
2002, CL2002 hereafter). 

In the space and luminosity distribution studies, it is very important to make corrections for the 
observational selection effects. In this contribution we used an empirical (and simple) method for the 
correction of observational selection effects, as in Kodiara (1974), Yusifov (1981), Leahy & Wu (1989), 
Wu & Leahy (1989, WL89 hereafter), with some modifications. 



2 Available Data and the Selection Effects 



At the time of preparing this contribution, the number of existing pulsars was 1412. The ATNF Pulsar 
Catalogue contains the data of PMPS, the Swinburne survey (Edwards et al. 2001) and all previous 
pulsar survey results, and it supplies a good sample of data for the statistical study of pulsars. Nearly 
600 of them were discovered in the Parkes and Swinburne Multibeam Pulsar Surveys (MBPS). 

In this study we are interested mainly in statistics of "normal" pulsars. For this reason we excluded 
from our sample binary and recycled (P < 10~^^ s/s), globular cluster and extragalactic pulsars. In 
our study we mainly used pulsars in the regions of the Galactic latitudes |&| < 15° and longitudes 
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— 100° < Z < 50°. In this study, as in the Galactic electron distribution model NE2001, we used the 
lAU recommended value Rq = 8.5 kpc as the distance Sun-GC. 

The apparent distribution of pulsars for the subsequent corrections due to selection effects is derived 
in the following manner. We drew equidistant concentric circles around the Sun (ri) and the GC (Rj) 
and made up a quasi-regular grid of points at the points of intersection of these circles on the Galactic 
plane. Then we drew circles of radius Rd around the grid points (or cells), counted the number of 
pulsars within the boundaries of circle Rd, and estimated surface densities of pulsars around grid 
points. Rci increasing linearly is chosen from O.li?0 to 1/10 of the largest considered distances from 
the Sun (18.7 kpc) and calculated by the relation Rd = 0.85 * (1.2 * (i — 1)/11 + 1), where i varies 
from 1 to 12. Similar variation of cell radii naturally reduces the density estimation errors at large 
distances due to distance uncertainties of pulsars. 

Making corrections for the selection effects, we used these data for the calculation of the radial 
distribution of pulsars in the Galaxy. Selection effects we divide into two categories and define them 
by the relation: 

p(r, R, l{r, R)) = K{l)K{r)po{r, R, l{r, R)) , (1) 

where I is the Galactic longitude; r and R distances from the Sun and the GC; K{1) direction and K{r) 
distance-dependent selection factors; po(r, R, l(r, R)) is the true and p{r, R,l{r, R)) is the observed 
distribution of surface densities of pulsars on the Galactic plane. 

K{1) is connected with the background radiation which leads to variations of the survey sensi- 
tivity with Galactic longitudes. Neglecting latitude dependence as in WL89, the direction-dependent 
correcting factor Ki{l) is defined by the relation: 

i^,(/) = i + Z^i-j|ziM, (2) 

where Tsky 1374(Z, is the sky background temperature at 1374 MHz and Tr = 21 K is the noise 
temperature of the Parkcs multibeam receiver (Manchester et al. 2001). The correction factor Ki(l) 
which was calculated from Eq. (2) is shown in Fig. 1. The observational number (or density) must be 
corrected for this selection effect, multiplying by Ki{l). It is evident that K{1) in Eq. (1) is equal to 
Ki^[l). In the anticenter direction and in the cell where the Sim is located, Ki{l) is assumed to be 1. 

The distance-dependent selection factor K{r) in Eq. (1) is the function of pulse broadening caused 
by scattering, scintillation etc. Here we make the simplifying assumption that the combined effect of 
these factors may be described by the exponential law as: 

K{r) = Ki{r)K2{r) ■ ■ ■ K,{r) = e.-"^'' e-"^'' ■ ■ ■ e""''' = e""'' , (3) 

where Ki{r), K2{r) etc. are the distance-dependent correction factors relating to the pulse broadening, 
scattering and other selection effects. The quantitative estimation of each of these factors indepen- 
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Fig. 1. Dependence of the correcting factor Ki{l) as a Fig. 2. Variation of the surface densities of pulsars on 
function of Galactic longitude. the galactocentric circle Rq = 8.5 kpc at various distances 

from the Sun. 
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dently is difficult. But the combined effect of these factors leading to reduced detection of pulsars 
away from the Sun may be estimated empirically as described below. 

Assuming that the surface density of pulsars is symmetric around the GC, and considering a 
galactocentric circle with the radius Rq, from Eq. (1) we obtain: 

p{r,RQ)=Kil)K{r)poir,RQ) . (4) 

where p(r, Rq) is the observable density; K{1) is known and Poir, Rq) is constant. 

In order to derive K{r) from Eq. (4), in Fig. 2, we have plotted p{r,RQ)/K{l) against the helio- 
centric distance by circles. 

Fitting the data in Fig. 2 by the least mean-squares (LMS) method for the distance-dependent 
correction factor we obtain: 

K{r) = exp (— cr) , (5) 
where c = 0.362 ± 0.017 and Inpo = 3.73 ± 0.16. 

3 Radial Distribution and the Luminosity Function of Pulsars 

The true densities of pulsars were found from the apparent densities, after correcting them by the 
relation (1). K{1) and K{r) are calculated from Eqs. (2) and (5) respectively. The radial distributions 
of densities from the GC may be found by averaging the densities p{r,R) at various distances {Rj) 
from the GC. 

The results of obtained radial distribution of surface densities of pulsars (squares) with correspond- 
ing error bars, together with the distribution of other Population I objects, are plotted in Fig. 3. In 
evaluating errors, we assumed that the number of pulsars in cells, in a rough approximation, follows 
the Poisson statistics. The errors for various data points are not equal, so we calculated the weighted 
average of surface densities and corresponding errors. In order to simplify the comparison with other 
results, the densities in Fig. 3 are normalized to the surface densities at the solar circle. 

Within the galactocentric radius 0.5 kpc pulsars are absent; and from the available data and 
precision of distance estimates, it is difficult to estimate the density of pulsars there. But within the 
range 0.5 kpc < R < 1 kpc there are 4 pulsars. Applying the correction factors K{1) and K{r) to the 
apparent density, we obtain 50 pulsars kpc~^. In fitting the radial distributions of pulsars, just this 
value will be used. 

The radial distribution of pulsar surface densities in Fig. 3 has been fitted by the relation: 

p(i?)=^(^Ayexp 

where Rq = 8.5 kpc is the Sun— GC distance. The best results of the fits are: A = 37.8±2.1 kpc~^, a = 
1.12 ± 0.10 and b = 3.20 ± 0.24. 

Revising Fig. 3, it is seen that high frequency, high sensitive PMPS searches revealed many pulsars 
around the GC, and within error limits the derived density (50±16 pulsars kpc~^) is not less than 
that in the circumsolar region. Applying the Tauris & Manchester (1998, TM98 hereafter) beaming 
model, for this value we obtain (500±150) pulsars kpc~^. 

In Fig. 3 we compare radial distributions of pulsars and SNR with the radial distributions of their 
progenitors. It is seen that, although the maxima of Population I objects coincide, the maximum of 
pulsar distribution is shifted to the GC by nearly 1.5 kpc. The maxima of pulsar and SNR distributions 
nearly coincide. But the radial scale length (RSL) of pulsar distributions (~ 4 kpc) is nearly two times 
less than the SNR distributions. 

In Fig. 3 we also plotted by dots the radial distribution of the birth location of NS which was derived 
by Paczynski (1990). Although the maximum of this distribution coincides with the maximum of the 
SNR distribution, its shape considerably deviates from the distribution of Population I objects. If the 
progenitors of NS are OB-type Population I stars, then from qualitative considerations of Fig. 3, it 
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Fig. 3. Radial distributions of pulsars (squares) and other types of Population I objects: SNR distributions are from 
Case & Bhattacharya (1998); H2 column densities from Bronfman et al. (1988) and Wouterloot et al. (1990); Hll regions 
are from Paladini et al. (2004) (it must be noted that Hll regions represent the number of sources in the 0.5 kpc wide 
Galactocentric rings, but not surface densities); OB star- formation regions data from Bronfman et al. (2000); marked 
by P are the radial distributions of birth location of NS (from Paczynski 1990); the expected radial distributions of 

birth location of NS (Eq. (7)) are marked by crosses. 



seems that the radial distribution of the birth location of NS must be located near the line shown by 
crosses in Fig. 3 and may be described by the relation: 



p(i?) = A(A) exp 



R. 



-b 



R 

Rq 



(7) 



where a — 4, b ~ 6.8 and Rq — 8.5 kpc is the Sun-GC distance. The constant A must be chosen from 
the calibration constraints (particularly, in the current case A = 1050). 

In order to construct LF, we selected the sample of pulsars within a cylinder of radius 2 kpc from 
the Sun. For faint pulsars we applied standard volume correcting methods, and LF is calculated as 

where is the true LF, ^o{L) is the apparent LF, V[< d{L)] is the volume of sphere with radius 

d{L) from the Sun in which pulsars with luminosities L are still observable, Vmax = 27rd^(L)zmax is the 
volume of cylinder of radius d(L) around the Sun and Zmax is the approximate value of the maximum 
height of pulsars from the Galactic plane. We accept that Zmax — 2 kpc, and pulsars outside these 
altitudes are not more than 2-5% of observable pulsars. 

The LF calculated on the basis of these data is presented in Fig. 4. The total number of observable 
pulsars in every luminosity bin is estimated on the basis of all available pulsars, but taking into account 
that high sensitive MBPS are carried out only in the region of |&| < 15° and —100° < I < 50°. 

For the low-luminosity pulsars the LF is estimated from relation (8), where Zmax is assumed to 
be 2 kpc. There are only two pulsars with luminosities > 100 mJy kpc^ in the considered region of d 
< 2 kpc around the Sun. <i>(L) values for them are calculated on the basis of large-distance pulsars 
revealed in the MBPS, assuming that they are observable almost everywhere in the Galaxy. 
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Fig. 4. Observed and corrected LF of pulsars. The observed number of pulsars in the luminosity bins are indicated under 
the corresponding bin. The observed distribution is given by the dotted line. A solid line shows the volume-corrected 
distribution of observable pulsars. The errors here represent the statistical uncertainty (30%) in distance estimates. A 
dashed line shows the distribution of active pulsars, after applying beaming correction. Lines marked by circles and 
crosses correspond to the Log— Normal fitting of observed and corrected data by the relation (9). A dash-dotted line 
corresponds to the fitting of the data above ~1 mJy kpc^ assuming a slope of —1. 



In this study we applied the TM98 beaming model for the correction of LF. We estimated approx- 
imate values of the beaming fraction (BF) corresponding to each luminosity bin and show them at the 
top of Fig. 4 in percent. LF, corrected according to these BFs, is shown by the dashed line in Fig. 4. 

The MBPS nearly doubled the number of pulsars around 0.1 mJy kpc^, but it is not enough to 
describe the LF by one straight line. As is seen in Fig. 4, only the LF in limited regions of luminosities 
may be fitted by a straight line. A flattening of LF at low luminosities was already noted in previous 
studies and explained as a result of the deficit of low-luminosity stars (see Lyne et al. 1985; Guseinov 
& Yusifov 1986; Cordes & Chernoff 1997; Lyne et al. 1998, and references therein). We fitted the 
data in Fig. 4 by the Log-Normal distribution before and after beaming corrections by one analytical 
relation: 

pulsars kpc~^ . (9) 

For the constants of distributions Al,(Tl and logLo before and after beaming corrections, we ob- 
tained values 67±15; 1.21±0.07; -1.64±0.33 and 2600±600; 1.24±0.10; -2.60±0.57 correspondingly. 



exp 



1 / logL - logLp 
'2 1 aL 



4 Results and Discussion 

We have studied the population of normal pulsars in the Galaxy, using the ATNF Pulsar Catalogue 
where the distances are calculated according to the new electron densities model NE2001. 

We have derived the radial distribution of surface density of normal pulsars by using the new 
distances. The maximum of radial distribution is located at ~3.5 kpc and ~1.5 kpc nearer to the 
GG relative to the maximum of distributions of Population I objects. Although the maximum of 
the distribution nearly coincides with the maximum of distributions of their progenitors (SNRs), the 
RSL of pulsars is ~4 kpc, i.e. nearly two times less than the SNR distribution. Integrating the radial 
distribution for the total number of normal pulsars with luminosities ii4oo > 0.1 mJy kpc^ in the 
Galaxy, we obtained (24±3) x 10'^ and (240±30) x 10'^ before and after applying beaming corrections. 
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The surface density of pulsars around the GC region is (50±16) pulsars kpc~^ and (500±150) 
pulsars kpc~^ before and after applying beaming corrections. 

For the luminosities at 1400 MHz, the LF of pulsars is constructed on the basis of pulsars within 
2 kpc from the Sun. In estimating beaming corrections we applied pulsar beaming models of TM98. 
The low luminosity end flattening of LF became increasingly more evident. The LF of pulsars may be 
better described by the Log-Normal distribution (Eq. 9). Above a luminosity value of 0.1 mJy kpc^ 
for the local surface density of normal pulsars we obtained (41 ± 5) pulsars kpc~^ and (520 ± 150) 
pulsars kpc~^ before and after beaming corrections. Assuming r = lO^'yr as the mean life time of 
active pulsars, for the corresponding birth-rate of pulsars in the solar neighborhood we obtain 4 and 
52 pulsars kpc~^ Myr^^ respectively before and after beaming corrections. 

We recommend a new relation (Eq. (7)) for the expected radial distribution of birth location of 
NSs, which is more closely related to the radial distribution of Population I objects. 

The total number and frequency of generation normal pulsars obtained from these data arc nearly 
1.5 and 4 times higher than previous estimates, correspondingly before and after beaming corrections. 
One of the reasons for the higher values is possibly connected with the overestimated electron densities 
of the ISM. This leads to the necessity of further improvement of the Galactic electron density model, 
which requires H i measurements for distant pulsars. 

One of the recent, but unsuccessful, attempts on pulsar distances and LF estimations are the 
publications of Guscinov ct al. (2002a, 2002b). However, due to the fallacy of their method and 
approaches to the problem and other numerous errors, their results are not discussed here. 

A more complete report on this work will be given elsewhere (Yusifov and Kiigiik, in preparation). 
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